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Abstract 

We formulate and study, in general terms, the problem of quantum system identification, 
i.e., the determination (or estimation) of unknown quantum channels through their action 
on suitably chosen input density operators. We also present a quantitative analysis of the 
worst-case performance of these schemes. 
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1 Introduction and background 

In quantum information theory all admissible devices are described mathematically by means 
of the so-called quantum operations (or quantum channels) [21 01 ■ 

Given a complex Hilbert space TC, denote by B(TC) the *-algebra of all bounded operators 
on 7i. In this paper we will work primarily with finite-dimensional Hilbert spaces, so that B(7i) 
includes all linear operators on H. Given Hilbert spaces H\ and Ti.2, a quantum channel T is a 
completely positive trace-preserving linear map of B(TCi) into 8(0.2) ■ All such maps admit the 
Kraus decomposition (Hj 



where : TLi — > TI2 are operators satisfying Ylk^k-^-k = ^-H\- This definition of the quantum 
channel is formulated in the Schrodinger picture, so that the density operators on 7i\ are mapped 
to density operators on Ti.2- The corresponding Heisenberg-picture definition goes the other way 
(observables on Ti.2 are mapped to observables on TLi) and yields a completely positive unit- 
preserving linear map T : B^H.2) B(7~ti) related to T by the duality 



for all p G BiTLi) and all X G 3(11.2). In this paper we will deal mostly with the Schrodinger 
picture. 

This seemingly simple framework turns out to be rich enough to cover all kinds of gen- 
eral transformations of quantum-mechanical states. In fact, both purely classical and hybrid 
(classical-quantum or quantum-classical) transformations can be included as well, simply by 
restricting to a suitable Abelian subalgebra either at the input or at the output. 

One of the basic challenges, both for theoreticians and for experimentalists, is to discover 
efficient procedures for analysis and synthesis of quantum channels. For instance, when designing 
a device for a specific task (e.g., an optimal quantum doner [1]), one has to run tests in order 
to determine whether the device performs according to specification. Several such procedures 
have been proposed already, such as the tomographic scheme of D'Ariano and Lo Presti [S] or 
the maximum-likelihood reconstruction method of Jezek, Fiurasek, and Hradil [6 . 

All of these schemes rely, in one way or another, one the one-to-one correspondence [7] 
between completely positive maps B{TL\) — > B^H.2) and positive operators on Ti.2 <8> Wi, to which 
we shall return later in this paper. Our purpose here is to phrase the ideas common to these 
schemes as an abstract problem of system identification. 

2 The quantum system identification problem 

Consider the following arrangement, shown in Fig.^ we are given a "black box" that implements 
an unknown quantum channel T : B(7ii) — > B(H.2), which we need to determine. This will be 
done by presenting to the black box certain suitably chosen input density operators p, thereby 
obtaining output density operators a = T(p), and then trying to determine (or to estimate) T 
given a set of ordered pairs (p,T(p)). 

We assume that we can re-use the black box any finite number of times, and that we can 
employ it as part of a more complicated arrangement. A typical strategy is to use quantum 
entanglement 0: one prepares input states (density operators) on the tensor product space 
T~ti ® Hi, and then subjects only one subsystem of the resulting composite system to T (see 




(1) 



A- 



Tr[T(p)X] = Tr[pT(X)] 



(2) 
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Figure 1: Input-output diagram for the quantum system identification problem. 
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Figure 2: The set-up for the quantum system identification problem that allows for the use of 
quantum entanglement. 

Fig.[2J). This way we have the channel T<g)id from BiTLi^TLi) into B(7^2<8>Hi), which corresponds 
uniquely to the original channel T. This extension is, in fact, at the basis of the Jamiolkowski 
isomorphism [2| (see next Section for a detailed discussion). 

The problem of quantum system identification can now be formulated as follows. Instead of 
viewing the arrangement shown in Fig.[2]as a mapping of density operators p to density operators 
T(g)id(p), we can think of it as a mapping of quantum channels T into density operators p[T] (we 
will use square brackets to distinguish the maps whose arguments are quantum channels from 
the maps that take density operators as arguments). That is, if we fix a density operator p, then 
we have p[T] := T <g> id(p). We will say that a density operator p is admissible if the map p[»] 
is invertible. Given an admissible density operator p, we will denote by p* the inverse mapping 
from density operators to channels. 

This points, at least in principle, toward a solution of the problem of quantum system iden- 
tification. All we need to do is to prepare an admissible state p, launch one of its subsystems 
through the black box T to get the output density operator a = T (g> id(p), and then reconstruct 
T as the inverse p^(cr). Of course now we are faced with (at least) two more problems. (1) 
What states are admissible? (2) What can we say about the performance of the reconstruction 
procedure as a function of the (admissible) input state? We will address these problems in the 
remainder of this paper. The rest of this section will be devoted to discussion of the general 
properties of the map />[•]. 

We would like to make some statements about the continuity of /?[•]. Let us equip the 
algebra B(TC) of linear operators on the Hilbert space Ti with the trace norm defined by 
||X||i := Tr|X|, where the absolute value of an operator X is defined as \X\ := \J X*X. Then 
\\A\\i = Tr A for any positive operator A, and furthermore \\p\\i = 1 for any density operator 
p. We will also need to estimate norm differences of quantum channels; the ideal norm for this 
purpose is the so-called norm of complete boundedness (or CB-norm, for short) JHJ) defined by 



(3) 
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where M n stands for the algebra ofnxn complex matrices. The norm || • || on the r.h.s. of (J3J) 
is the operator norm, defined for a general linear map M : B{H\) — > B^H.2) by 

||M|| := sup \\M{X)\\ X . (4) 

XeB(Wi);||X||i<l 

(Note that the above definition is tailored specifically for quantum channels in the Schrddinger 
picture; consult the monograph of Paulsen jlUj for generalities.) We have ||T|| c b = 1 for any 
quantum channel T . We shall have an occasion to use some other properties of the CB-norm in 
later sections; all we need right now is the inequality ||T(A)||i < ||T|| c b||A||i, which is obvious 
from definitions, and the multiplicativity of the CB-norm with respect to the tensor product, 
US' <8> T|| c b = HSHcbllTllcb. 

With these preliminaries out of the way, consider a density operator p. We can easily see that 
the map p[»] is continuous. Indeed, consider two quantum channels, T and T 1 . By definition we 
have 

Wpm-pn = nr®id(p)-T'®id(p)ii 1 

< ||T(g>id-T / <g)id|| c b 

= ll^-T'Heb (5) 
(we have used the fact that the trace norm of a density operator is equal to one). 



3 Admissible states and the Jamiolkowski isomorphism 

In this section we describe an approach to the construction of admissible states. This method 
is closely connected to the Jamiolkowski isomorphism 7 j, and is, in fact, a natural extension of 
the latter. First we need some mathematical machinery from the theory of completely positive 
maps. 

Given a C*-algebra ^4, denote by A + the cone of positive elements of A (i.e., precisely those 
elements that can be written in the form A* A for some A £ A). A linear map T between C*- 
algebras A and B is called positive if T(A + ) C B + , and completely positive (CP, for short) if the 
maps T <g> id n : A (g> M n -> B <g> M n are positive for all n £ N. A quantum channel is thus a 
specific instance of a CP map. 

According to the fundamental theorem of Stinespring for any CP map T : B(H\) — > 
B(TL2) there exist a Hilbert space £ and a bounded operator V : H2 — > T~ti <S> £ , such that 

T(p) = V*{p®l £ )V. (6) 

The pair (V, £) is called the Stinespring dilation of T. Furthermore, with the additional require- 
ment that the linear span of the set {(A <g) lg )Vij) \ A G B(Tti),ip € TL2} be dense in Hi £, the 
pair (V, £) determines the CP map T uniquely up to unitary equivalence. In that case we speak 
of the minimal Stinespring dilation of T. 

Consider now the set of all CP maps between B{TL\) and B(H.2), for some Hilbert spaces Hi 
and H2- This set can be partially ordered in the following way. Given two CP maps T\ and T2, 
we will say that T\ is completely dominated by T2 JT21 (and write T\ < T2) if T2 — T\ is also a CP 
map. If T\ < XT2 for some positive A £ R, we will say that T\ is completely \-dominated by T2. 

Given a CP map T, the set of all CP maps completely dominated by T can be characterized 
completely using a theorem of Belavkin and Staszewski , which is referred to as the "Radon- 
Nikodym theorem" for CP maps, and asserts the following. Let (V, £) be the minimal Stinespring 
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dilation of a CP map T. Then a CP map S is completely A-dominated by T if and only if it has 
the form 

S(p) = V*(p®F)V (7) 

for some positive operator F E B{£) with \\F\\ < A. Furthermore, the operator F is determined 
uniquely by S and (V,£). 

Next we would like to show that the "Jamiolkowski isomorphism" between CP maps of 
B(Hi) into B(H.2) and the positive operators on H2 <8> Hi is a direct consequence of the above 
Radon-Nikodym theorem |13| . 

We will consider quantum channels from S(Hi) into B(H.2), where Hi and H2 are finite- 
dimensional Hilbert spaces. Let us fix an invertible density operator p on Hi, which we will call 
the reference state. Let and be the eigenvalues and the eigenvectors of p, and let us also 
fix an orthonormal basis {f p } for H.2- Denoting by £ the tensor product H2 <8> Hi, define the 
isometry Vp : H2 — ► Hi (8> £ by 

V p iP:=Y,p] /2 {fM)<i>i®fn®<l>i- (8) 

1/2 

Consider now a channel T. Let us define the unit vector Q p E Hi <8> Hi by Sl p := YliPi 4*% ® & 
and the positive operator on £ by 

F TjP := (1 ® p- x )T ® id(|« p ><n p |)(l /T 1 ). (9) 

With these definitions, we can write 

r(«7) = F;(<7®F r , p )^, (10) 

where the action of the coisometry V* on the elementary tensors £ (8) € Hi £ is given by 

®v) = J2p¥ 2 (U ® ^><<&IO/„, (11) 

t./i 

and then extended to all of Hi £ by linearity. 

It is now an easy consequence of the Belavkin-Staszewski theorem |12| that the operator Ft p 
uniquely determines the channel T, and that for any positive operator F G B{£), the map 

M{a) = v;(a®F)V p (12) 

is completely positive. We see therefore that any invertible density operator p on Hi gives rise to 
an admissible pure state to = \Q p ){Q p \, in the sense that the mapping u[T] := T®id(\£l p )(tt p \) is 
invertible. That is, the image of any density operator w on £ under the inverse map p* is given 
by 

pKw) = V* (• ® (1 ® p~ x )w{l ® p- 1 )) V p . (13) 

It is important to realize that, in general, p\w) is not a quantum channel, unless w satisfies the 
additional consistency condition 

Tr % [(1 ® p-^Ml ® P' l,2 )\ = l«a • (14) 

This will hold automatically in the quantum system identification setting (see Fig. |2J) , provided 
that there is no additional noise in the apparatus. We note also that the Jamiolkowski isomor- 
phism is a special case of this formalism ^3]> and is obtained if we pick as the reference state 
the maximally chaotic density operator (dimHi) _1 l-^ 1 . 
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Before we go on, let us remark that the set of pure states on Hi<8>Wi obtained by "purification" 
of invertible density operators on 7i\ does not exhaust all possibilities for admissible states. In 
a recent paper, D'Ariano and Lo Presti ^1] have constructed a wide class of admissible states, 
which includes subset the states discussed here. 



4 The performance of quantum system identification procedures 

In this section we will quantify the performance of quantum system identification procedures as 
a function of the admissible state used as an input to the unknown channel. 

Consider two density operators w\ and W2 on H2 <8> *Hi that satisfy the consistency condition 
(|14jl. Then there exist quantum channels T±,T2 : B(Tti) — > B{TL2) such that 

w i = p[T i \=T i ®id(\n p )(n p \), i = l,2. (15) 

Furthermore, from (|13|) it follows that 

T i = P *(w i ) = V;(.®F i )V p , (16) 

where := (1 (g> p~ 1 )wi(l (g) p -1 ), « = 1,2. To determine how close the reconstructed channels 
p^(w\) and p"(u>2) will be when the corresponding density operators w\ and W2 are close (say, in 
trace norm), we will get a lower bound on the channel fidelity ^3] between pfi(wi) and p*(u>2)> 
defined in the following way. Consider two channels T\,T2 ■ B(7ii) — > B(Tt2), and define the 
density operators 

ai:=Ti®id(\n){n\), (17) 

where Q := (dim^i) -1 / 2 £^ (g) e^, the summation taken over some orthonormal basis of 7i\. 
Note that |^)(^| is an admissible state corresponding to the maximally chaotic density operator 
(dimWi)" 1 !^ . Then the channel fidelity ^Hl is defined by 



HTl,T 2 ):= (TrV^WPj , (18) 

where the quantity on the r.h.s. of ()18|) is the mixed-state fidelity |16j . We do not need all of the 
properties of the channel fidelity ((TH|) [but see Ref. HS|], except the following: 

2-2(^(T 1 ,T 2 )) 1/2 < Iki-^Hx, (19) 

which is a simple corollary of the results of Fuchs and van de Graaf |17j . We also note that the 
channel fidelity has the natural property that J-(Ti,T2) = 1 if and only if T\ = T2 (this is a 
straightforward consequence of the properties of the mixed-state fidelity |16j). 
We can rewrite w\ and W2 from (|15|) in terms of p, ai, and <r 2 ■ 

(1^-V 2 H(18 P -V 2 ) . 
* = dkn^ ' , = 1 ' 2 ' (20) 

Then we can use the well-known inequalities ||AB||i < and ||^4*A|| = ||A|| 2 , where || • || 

is the usual operator norm [5], to get 

hi-*2\\i< • (21) 
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Combining this estimate with (|19[). we obtain 



2-2y/r( f fl(w 1 ),ffl(w 2 )) 

I ' 11^1 ' 
dim Tii 



Upon rearranging, we get the desired lower bound: 

T(p i {wi),p i (w 2 )) 



2 



K i -2iyk |K - , " 2iii j ■ w 

We see right away that worst-case performance of the channel reconstruction procedure is con- 
trolled by the smallest eigenvalue of p (or, equivalently, by the largest eigenvalue of p~ l ). This 
fact has also been pointed out by D'Ariano and Lo Presti [Sj, and Eq. Q23JI gives the correspond- 
ing quantitative estimate. Note that in the case of p = (dim "Hi) 1?^, the constant in front of 
the trace norm on the r.h.s. of Q23J) is 1/2, which yields worst-case performance that depends 
only on the states w\ and w 2 - 

Note that we have discussed here the ideal scenario, namely that there is no additional noise 
in the apparatus used for the channel reconstruction. Any such disturbance will, of course, 
further degrade the performance of the scheme. 



5 Discussion and conclusions 



We have outlined a general mathematical framework for quantum system identification, i.e., the 
determination (or estimation) of quantum channels through their action on suitably chosen input 
density operators (we have called them admissible states). In general, the channel reconstruction 
procedure will involve the preparation of an entangled state, followed by the application of an 
unknown channel to one of the subsystems, leaving the other one intact. One can show [51 1141 US] 
that the use of entangled states results in an overall improvement, in either the precision or 
the stability of the reconstruction procedure. On a more fundamental level, however, the use 
of entanglement is also essential in light of the one-to-one correspondence between quantum 
channels and bipartite density operators [that satisfy the consistency condition ([140]. which can 
be explained in abstract terms within the framework of Radon-Nikodym type theorems for CP 
maps PIES]. 

In this paper we have emphasized quantum channels acting on finite-dimensional algebras, in 
order to keep the presentation simple. However, it is important (e.g., for quantum information- 
theoretic studies in quantum optics) to have a mathematical theory of quantum system identi- 
fication in infinite dimensions. Some steps in this direction have already been taken (see, e.g., 
D'Ariano and Lo Presti ^1] or Raginsky J3]). Let us briefly comment on some of the big points. 
One starts with a density operator p that is invertible; however, the inverse is now an unbounded 
operator. This implies that the reconstruction map will fail to be continuous, which will result 
in an unbounded growth of statistical errors during the tomographic estimation of the matrix 
elements of the Radon-Nikodym density Ft, p - 
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